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Abstract 



We propose a phenomcnological three component model describing the decay amplitude 
of the process Df — > 3tt. The first component is a constant background, the second one is a 
Breit-Wigner amplitude associated to a quasi two-body /o(980)7r + state and the third one another 
Breit-Wigner amplitude corresponding to a possible quasi two-body p{71Q)iT state. We show that 
it is possible to reproduce the observed total rate for Df — > 7r + 7r + 7r _ as well as the two other 
measured branching ratios for the non resonant part and the resonant foir + one. Implication of 
the pir + state, of which an experimental limit has been given, is discussed. An upper bound in the 
10 MeV range for the decay constant f n > of the 7r(1300) meson is obtained . 

Predictions are given for the Df — ► tt tt tt + rate as well as for the ir + and n~ (ir° and 7r + ) 
energy distributions for these two decay modes Df — ► 7r + 7r + 7r~ (£)+ — > 7r°7r°7r + ) respectively. 



PACS numbers : 12.15.Ji, 13.25.Ft, 14.40. Aq, 14.40.Lb 
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I Introduction 



We are interested, in this paper, in the decay of the Df meson into three pions. Experimental 
data are available only for the final state ir + ir + ir~ ||, |j, but not for the 7r°7r°7r + one. By inspection of 
the results, we observe that the quasi two body state /o(980)-7r + has been clearly identified [0 whereas 
for the quasi two body /9°(770)7r + state, no significant signal has been found Q. Another remarkable 
feature of the data is an important fraction of the decay proceeds into a non resonant 7r + 7r + 7r _ state 
[0, this fact necessarily implies a sizable W + annihilation contribution. More precisely, the rate for 
.Dj" decays into a non resonant 7r + 7r + 7r~ state is about 1/3 of the dominant spectator decay rate 
— > (j)ir + , both being governed by the same Cabibbo-Kobayashi-Maskawa (CKM) favored factor 

V cs V ud- 

After a brief discussion of the kinematics for three body decay and a presentation of the experi- 
mental results in section II, we construct a three component model for the transition amplitude of the 
mode Df —* 7r + 7r + 7r~. The first component corresponds to a quasi two body state, followed by 
p° — > 7r + 7r~. Since p° meson is a (uu — dd)/\^2 isovector state, the decay Df — > p°7r + can proceed 
only through the W annihilation mechanism (henceforth WA). In section III, we propose a model for 
the W + — > p°ir + transition involving an intermediate virtual 7r(1300) meson. The experimental upper 
limit for the branching ratio B(Df — > 

A + 

) allows us to derive an upper bound for the decay constant 
f i of the 7r(1300) meson to be around 10 MeV. 

The second component is associated to the observed quasi two body /ovr + state followed by the 
decay /o — > 7r + 7r~. The quark content of the /o(980) meson is crucial for constructing the amplitude 
.D+ — > /ovr + . If the fg is a pure (uu + dd)/\/2 isoscalar state, we must again use the WA, as for the 
p° previously. In section IV we show that by comparing the — > p°7r + and — > /ovr + rates, such 
an assig nment (uu + dd)/y/2 for / (980) is in strong disagreement with experiment. Then the /o(980) 
must content a large ss component or even be a pure ss state. Now the spectator mechanism can be 
used for the decay Df — > fon + , and from experiment, we estimate in section IV the magnitude of the 
hadronic form factor Fq s ^° (m^) describing the D s — > foW + transition. If /o(980) is a pure ss state, 
we obtain, from experiment i ? ( ^ >s ^ (m^) ~ 0.36, a reasonable value. 

For the third component associated to a non resonant 7r + 7r + 7r~ state, we introduce a complex 
constant amplitude which corresponds to an uniformly populated Dalitz plot. Such a component is 
determined, from experiment, in section V where, in addition, the total rate for Df — > 7r + 7r + 7r~ is 
computed using the three component model described above. The result, consistent with experimental 
data, is presented as a function of a possible phase for the constant non resonant amplitude. 
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The 7r + and it~ energy distributions for the decay Df — > 7r + 7r + 7r~ are computed in section VI. 
The /ott + quasi two body state is clearly seen in the ir + energy distribution, as already obtained in 
|ffl], no p°tt + enhancement is present, and the tt~ energy distribution does not reveal any particular 
structure even if it is very different from the pure phase space distribution. It is then straightforward 
to make predictions for the decay mode Df — ► 7r°7r°7r + using isospin considerations relating the two 
final states 7r + 7r + 7r~ and 7r°7r°7r + . In section VII, we present our results for the rate and for the 7r°, 
7r + energy distributions of the decay into 7r°7r vr+. 

Finally, we end this introduction by few considerations concerning the quark content of the /o(980) 
scalar meson based firstly on SU(3) flavor symmetry and secondly on the decay products of the /o(980). 
The SU(3) flavor nonet J PC = ++ contents an isovector meson ao(980), two isoscalar mesons /o(980) 
and /o(1300) and two isodoublets of strange mesons 1^(1490). The physical particles /o(980) and 
/o(1300) are expected, as usual, to be the result of a configuration mixing between the isoscalar octet 
and singlet weights. We remark that the quasi degeneracy in mass between the ao(980) and /o(980) is 
reminiscent of that between the /o(770) and o;(782) in the J PC = 1 nonet. A first naive expectation 
is that / (980), like w(782), is dominantly a (uu + dl)/y/2 bound state, the / (1300) like <£(1020) 
being essentially a ss state. However such an assignment is wrong for at least two reasons from both 
theoretical and experimental considerations. Indeed if one tries to apply for the scalar J PC = ++ 
nonet the usual mixing formalism which has been very successful in the J PC = 1 , h , 2 ++ , 3 
cases, then we are faced with the impossibility to determine the mixing angle 9 S for scalar mesons. 
Because of the large mass of the 1^(1490) meson, the formalism leads to tan 2 9 s < 0, no matter 
is used the quadratic or linear Gell-Mann - Okubo mass formula. Therefore the previous analogy 
between the J PC = 1 and ++ nonets is meaningless. The second reason is an examination of 
the experimental decay products of the /o(980) and /o(1300), in particular the -kit and KK channels. 
In spite of the tiny phase space available for /o(980) — > KK, its branching ratio 22% is actually a 
very large number. The situation is inverted for the second meson /o(1300), the available phase space 
for KK is much larger, however its branching ratio is only 7.5%, with 93.6% into tttt. From these 
experimental considerations, we conclude that if the J PC = ++ nonet is made of quark-antiquark 
bound states, the /o(980) must contain a large ss component and the /o(1300) a small ss one, with a 
(uu + dd)/y2 component being dominant for /o(1300) and secondary for /o(980). Moreover, we shall 
see in section IV that the comparison of the D s — ► p ir + and D s — ► /o(980)7r + rates also confirms the 
dominant ss component of the /o(980) meson. Such ss classification for /q(980) has also been recently 
advocated in ||. 
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II Generalities and Experimental results. 
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1). We study the decay of the Df meson of energy momentum P into three ir mesons of energy 
momenta p±, p 2 , P3 with the relation P = pi + P2 + P3- We introduce the Mandelstam variables 
si,s 2 ,s 3 and the ir meson energies in the Df rest frame E\, E 2 , E%. At the level of accuracy of the 
experiments, it is legitimate to neglect the mass difference between charged and neutral pions and we 
get 



(1) 



Energy momentum conservation implies the relations 

si + s 2 + s 3 = m 2 Ds + 3ml, E 1 + E 2 + E 3 = m Da (2) 

The double differential distribution is given in terms of the transition matrix element by the 
expression 

dT = — ^ — | < 37r|T|D+ > | 2 dE x dE 2 (3) 
647r d mo, 

We remind that the transition matrix element < 3ir\T\Df > involving three body final state is 
dimensionless. In the (Ei,E 2 ) plane, the phase space is defined by the constraints 

(m 2 D - 3ml) 
2m n a 

E-(Ei) < E 2 < E+(Ei) (4) 



with 



E±(E) = l (m Ds -E)± \ 



(E 2 — m 2 T )(m 2 Ds — 3m 2 — 2mD s E) 



(5) 



2 V Us ' 2\ ml s + ml-2m Ds E 

Of course, the mass difference between charged and neutral pions being neglected, we have the same 
phase space in the two other planes (E 2 ,E 3 ) and (Ei,E 3 ). Two possible three pions final states are 
7r + 7r + 7r~ and 7r°7r°7r + . We make the following choice of ir meson variables 

(1) . final state tt+tt+tt-, E 1 (ir + ), E 2 (ir + ), E 3 (ir-) 

(2) . final state 7r°^°^+, E^n ), E 2 (tt°), E 3 (tt + ) 

2). Let us now discuss the experimental situation concerning the only observed 3 pion final 
state 7r + 7r + 7r~. The decay rate for the mode Df — ► 7r + 7r + 7r~ has been first measured by the E 691 
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collaboration at Fermi Lab. [Q| and the result has been confirmed few years later within error by 
the NA 82 collaboration at CERN In addition, the E691 group has been able to exhibit in the 
7r + 7r + 7r~ Dalitz plot a structure due to the fo(980)n + quasi two body state. However a structure due 
to the /9 (770)7r+ quasi two body state has not been clearly identified. The various branching ratios 
in are normalized to the most easily measured mode Df — > <pir + . Using the Particle Data Group 
value Q 

B(D+ -^(f> + ir + ) = (3.5 ± 0.4)% (6) 
the branching ratios obtained in Q are : 

B(D+ -> n + TT + TT-) TO T = (1.54 ±0.42)% (7) 

B(D+ ->7r + 7r + 7r~)iVfl = (1.015 ± 0.352)% (8) 

B(D+ -► / vr + ) = (0.98 ±0.38)% (9) 

B(D+ -> p°vr+) < 0.28% at 90 % C.L. (10) 

where, in (^) the subscript A^i? means non resonant Wi. Our results (@), (S), and (|l^) agree with those 
given in Q. However the branching ratio of Eq.(^) quoted in Q is slightly smaller, (1.35 ± 0.31)%, 
because an average has been made between the data obtained in Q and @. For a question of 
consistency we shall use, in this paper, the set of data (]?]) — ( |l0| ) in order to preserve the relative 
fractions of the various terms as quoted in 
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III. The decay mode for Df — > p°n + — > tv + tt + tt 

1) . Let us consider the quasi two body final state p°7r. The transition matrix element has the 
form 

< A + I T | D+> = g Ds(m (P Ds +PnV e*(p p ) (11) 

where e p (p p ) describes the p meson polarization. The dimensionless quantity gD 3 pn is a weak complex 
coupling constant in terms of which, the decay width is given by 

1 K 3 

T(Dt^p°n + ) = - ^ \g Dsfm f. (12) 

Here K p is the center of mass momentum of the final pn particles in the meson rest frame. 

2) . The p° and pion are quark-antiquark bound states containing only u and d quarks and the 
decay Df — > p°ir + can only proceed through the WA, since both charm and strange quarks being 
absent in the final state. For the transition W + — > p°7r + we assume, by the Partial Conservation of 
the Axial Current (PCAC) the contribution of an intermediate state having the quantum number of 
a 7r + meson M, ^, |?j and which might be the tt meson itself or its recurrence 7r(1300). It has been 
checked that the n meson intermediate state gives a contribution many order of magnitude smaller 
than that of the 7r(1300) and only the latter one is retained in the diagram of Fig.l. 

Here tt' = 7r(1300) and the corresponding expression of the weak decay constant gD a p-K is given by 

Gf m 2 D # f Ds f i 

9d sP -k = «l /=— 1 V cs v ud — 2 n "p ( — 2~T 9-k'ptv ( 13 ) 

V2 mr, — mi, -imd ' mi ^ 

where ai is the phenomenological parameter introduced by Bauer, Stech, Wirbel (henceforth BSW) 
Hi) the relevant CKM matrix elements; fjj a and jfy are the leptonic decay constants of the 

and tt' mesons respectively. 

The dimensionless strong decay constant g n ' is related to the corresponding decay rate by 

where K p is the C. M. momentum in the ir rest frame of the pn final state. 
Combining Eqs. ©, (0) and @ we obtain 

3 

(15) 



B(D+-A + ) = ^ (fj) # A 2 - A + ) 



where A is given by 



3). For masses, widths and CKM matrix elements we use the values collected by the Particle Data 
Group M. The quantity fjj a is choosen to be fjj a = 280 MeV consistent with recent experimental 
data || and theoretical expectations. For the decay 7r(1300) — * pn we use j|] 



B(tt — ► pn) 



0.3205 



(17) 



and we leave / / as a free parameter. For a% we take a\ = 1.26 ||. 

In the following, we use the simple energy dependence mT(s) = y/s T(m 2 ) for the widths, although 
more sophisticated expressions have been proposed [0, ||]. Because of the not too large difference 
between the 7T and masses, we expect the sensitivity of Eqs.(13) and (|l6| ) to different forms of 
T^i(m 2 D: _) to be relatively modest. 



We now write the branching ratio B(D S 

B(D+ 



p°ir + ) in the form 
= C 2 f, 

V TV 



(18) 



We retain only in the numerical computation of C 2 the large uncertainty due to the 7r(1300) width 
(J, T / = (400 ± 200) MeV and we get 



C' 1 



16.3 t \\ GeV~ 



(19) 



Using now the 90% confidence level upper limit B(D S -> p°vr+) < 0.0028 of Eq.(0) we deduce an 
upper limit for / / 



4' 



< 



11.45 
13.11 
17.72 



MeV for T < 

J 7T 



600 
400 
200 



MeV 



(20) 



Theoretical value for f w i, as large as 40 MeV, seems to be overestimated in the literature M. Eq. (pC|) 
gives, to our knowledge, the first experimental information on f t, and our result is compatible with 



theoretical calculations using non-relativistic chiral quark model [ 10 1 or QCD sum rule technique [11|. 

4). We now take into account the non-zero p meson width effect. The production of the quasi two 
body state p°ir + is followed by the decay p° — ► ir + iT~ . For the final state it + it~it + , E\ and E2 are the 
7r + energies and E3 is the tt~ energy. The decay amplitude for Df — > p°7r + — > it + -k~it + is given by 

' si - s 3 I 



T I Df > 



9D B pir 9pirir 



S2 ~ S3 



+ 



m 2 -si-im p Tp(si) 
and the corresponding decay rate has the form 



r(Z>+ /9°7T + ^ 7T+7T 7T+) = m Ds \g Da pn\ 2 I fipTTTr | 2 ^ 



s 2 -im p T p (s 2 



(21) 



(22) 



where J is the phase space integral 
1 



J 



m 



D 3 



dEi dE 2 



S2 ~ S3 



mi 



si - i y/s~[T f 



+ 



si - s 3 



111 



S2 



S2 r. 



(23) 



The energy dependence of the p meson width has been taken as explained previously and the factor 
l 

2 



1 in Eq.p^) is due to the presence of two identical tt + in the final state. 



The dimensionless constant g p7rn is related to the p meson width by 

iV 



7T 7T 



c o \ypiT1T ) 

D7T m z p 



(24) 



where is the C. M. momentum in the p° meson rest frame of the two pions. 
Combining Eq.(B^) and (1241) we obtain 



IT IT TT 



3 m Ds 2 T(p° -» 7T+7T ) 

\9D sP n\ (1 _ 4m 2 /m 2)3/2 



J 



ir* rrir, 



(25) 



It is straightforward to estimate the importance of the non-zero p meson width. The result is obviously 
independent of gD a pn and comparing Eq.([l2|) and (|22|) we get 

T(Df — ► p°ir + — > 7r + 7r~7r + ) 3 Ta p rar) a T(p° — > 7r + 7r~ 



T(Df -> p°7r+) 4vr ^3 (1 - 4m2/m2)3/2 

The value of J obtained by numerical integration is J = 8.0997 and we get 

F(D+ -> / o°vr+ -» -K+1T--K+) 



J 



T(DS 



p°ir^ 



0.9719 ±0.0096 



(26) 



(27) 



where the error is due to the uncertainty of the p° meson width. 

Finite p° meson width effect is small and the upper limits obtained in Eq. (20) for |/ /| are only 
slightly modified 



< 



11.61 
13.30 
17.98 



MeV for T , 

J 7T 



600 
400 
200 



MeV 



(28) 



IV. The decay mode for Df — > f ir + — > 7r + 7r + 7r 

1). We consider now the quasi two body final state /o(980)-7r + . Since two body decay amplitude 
has a mass dimension, let us write the transition matrix element for spinless particles in the form, 



< f ir + | T | D+ > 



m Ds 9D s f n 



(29) 
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where the dimensionless quantity gD s f -w is a weak complex constant in terms of which, the decay 
width is given by 

r( D+ - /o tt+) = ^K f \g Dsf ^\ 2 . (30) 

Here Kf is the C. M. momentum of the final particles /ovr + in the rest frame. 

2). The considerations made in the introduction concerning the quark content of the /o(980) 
meson can be tested in the decay of interest here, Df — * foTT + . Let us first assume that /o(980) is a 
(uu+dd) I y/2 bound state, then we can show that this assumption is in contradiction with experimental 
data. Indeed under such assumption, the only possible mechanism for both decays D s — ► fon + and 
D s —* p°7r + is the WA and the coupling constant gD a f n is obtained from a diagram similar to the one 
of Fig.l where the p° is replaced by the /o- It is straightforward to derive a relation between the two 
weak coupling constants 

where the subscript WA indicates that the decay amplitude is performed in the W annihilation model 
of Fig.l, from which we get 



T(D+ -> f 7r+) WA (K' p y (K f \ ( m' n \ 2 r(vr / f ^ 



(32) 



r(£>+^p°7r+) \K p ) \K' f ) \m Ds ) r(7r'^p°7r- 

where K'^ is the C. M. momentum in the ir' rest frame associated to the decay ir — > /ott + . 

We assume that the experimental result of the tt decay given in Q for the final state (vr7r)s_„, a „ e + 
7r + corresponds to fon + and we use 

r(7r' / vr+) 



r(7r' — » p°ir^ 
Numerically, we obtain from Eq.(|3z 



2.12. (33) 



T(D+ -» f 7T + ) WA = 0.318 • F(D+ -» p vr + ). (34) 

Using the upper limit ( |l0|) for the rate of Df — > p°ir + , we get the prediction 

B(D+ / 7T+)^ < 8.9 • lO" 4 (35) 

in complete disagreement with the observed experimental value Eq.(||). 

It follows that the dominant component of /o(980) cannot be (uu + dd)/y/2 and a large ss compo- 
nent must be present. In fact such an ss component will be responsible for the observed strong decay 
mode fo — ► KK. 

3). Having shown that /o(980) is mainly a pure ss state by both considerations based firstly on its 
decay modes into KK, secondly on those of Df into /o7r (beside the mixing in the J PC = ++ SU (3) 



flavor nonet argument mentioned previously in the Introduction), the dominant decay mechanism for 
the mode Df — > /ovr + is then the spectator diagram. The weak coupling constant gD a f n is now 
computed from the diagram of Fig. 2 and the result is 

G F m 2 Ds {m 2 Ds - m 2 fo ) f n Dafo 2 

ODsfoK - i a\ j=— Vcs Kd —2 — — ( m J \ 6b ) 

Vz m Ds m Ds 

where F Ds ^°(q 2 ) is the D s —* fo hadronic form factor in the BSW notation associated to the spin zero 
part of the weak axial vector current. Using the experimental branching ratio B(Df — > fo^ + ) in Eq. 



), we deduce the value of the hadronic form factor F^* * 



2> 



F^{mi) = 0.36 1 (37) 

which appears to be an acceptable order of magnitude for an hadronic form factor. 

If the /o(980) is not a pure ss state, what has been obtained numerically in Eq.(37) is the product 



of Fq s ^(m^) by a mixing coefficient /i corresponding to the amount of ss contained in /o(980). Of 
course, |/u| < 1, such that the true form factor is \F Ds ^° (m 2 ,)] > 0.361 oog- 

It is interesting to compare this phenomenological value with the similar one coming from the well 
measured Df —* 07r + decay mode. In the ideal mixing situation close to reality, the (ft meson is a pure 
ss state and the decay Df — > (f)TT + is described by the spectator mechanism involving the hadronic 
form factor Aq s< ^ '(m 2 ,). Using the experimental result quoted in Eq.(^) we obtain an experimental 
value for A^ a ^{m 2 ) 

\A° s<t, (ml)\ = 0.72 ± 0.04 (38) 

Of course this form factor is larger than the one previously obtained for ji F Ds ^° (m 2 ), but there are 
no reasons why they should be equal. 

However let us remind that both F Ds ^° and Aq s< ^ are determined by the factorization method. It 
is possible that non factorization contributions will change these numerical estimates in a sizable way, 
unfortunately, to our knowledge, a convincing approach to compute the non factorization term is not 
yet available. 

4). We now take into account the fo meson width effect. The decay amplitude for — > /o 7r+ — * 
ir + ir~ir + is obtained with the help of the definition Eq . (p9|) , together with the similar expression 
associated to the decay fo — ► ir + ir~. 

< IT 7T 7T | 1 | U s >/ 07r = gD s f 7T 9f TT+n- 



m% m\ 



m 



- si- i m fo T fo (si) ml - s 2 - i m fo T fo (s 2 ) 



+ —2 1 T^—F^TT (39) 
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and the corresponding decay rate is 

T(Df - / 7T + - 7T+7r-7T+) 



1 



if 



m l> s I |2 2 _ 

64 7T 3 m 2 l9Dsfonl 9 f^~ 2 

JO 



(40) 



where the dimensionless coefficient K is the phase space integral 



K 



m 



m 



dEi dE 2 



+ 



in 



fa 



si - i yfsl T fo ' m 2 fo - s 2 -i ^ F fo 



(41) 



The energy dependence of the /o meson width has been taken as explained previously and the 
factor | in Eq. ([!(]) is due to the presence of two tt + in the final state. 

The dimensionless constant gf 07T + n ~ is related to the partial decay width by 



r(/o -» tt+tt ) = — K w ( g fon+7T - f 

where K n is the CM. momentum of the two pions in the /o rest frame. 
Combining Eqs.(]40|) and fl42] ) we obtain 

f \ 3 



(42) 



T(D+ - / 7T + - 7r+7r-vr+) 



1 



8^ 2 



r (/o -> tt+tt ; 

1 — 4m 2 / m 2 



# lto s / 7r| 



(43) 



.fo 



The effect of the non zero /o width is estimated by comparing Eqs.(|30|) and 

\ 3 „ 



1 



fo 



K 



l-4ml/m 2 fo K f 



(44) 



T{Dt - / 7T+) X B(/ - 7T+7T-) 7T ^ m /o ^ 

The value of IT obtained by numerical integration is X = 5.7217. We get 

r(£>+ f 7r+ -> 7r+vr-7r+) . . 

- \ s ^- '— = 0.9886 ±0.1893 (45) 

T(D+ -> / 7T+) X B(/ - 7T+7T-) V ' 

where the error is due to the uncertainty on the /o meson width taken to be Tf Q = 47 ± 9 MeV. 

The /o finite width effect is small and the value of F^ s ^° (m 2 ^,) obtained in Eq.(p7|) is only slightly 
modifed. 



V. Total rate for the decay D+ — > tt+tv+tv 

1). We propose a model where the decay amplitude < it + it~it + \ T \ Df > is written as a sum of 
three contributions 

< 7r + vr"vr + \ T \ D+ > = < 7r + 7r"7r + | T \ D+ > p7T + < 7r + vr"vr + | T \ D+ > fo7T 

+ < 7r + vr"7r + | T | £>+ >nr (46) 
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The two first components in Eq.([46|) are associated to the quasi two body final states p°7r + and fo^ + 
and their expressions have been given in Eqs.(21) and (|39|). The third component describes a non 
resonant 7r + 7r _ 7r + state and is naturally assumed to proceed through the W + annihilation mechanism. 
In this case the general structure of the decay matrix element is 

< 7r+7r-vr+ \T\D+ > WA = i ai ^ V* V ud f Ds < | A„ | > (47) 

The matrix element of the divergence of the weak axial vector current between the vacuum and the 
three pion final state is an unknown structure function F(Ei,E 2 ) depending on two independent 
variables, choosen as the pion energies. 

Pg s < 7r+7r-7r + | A» | > = m Ds F(E 1 ,E 2 ) (48) 

Since the matrix elememt of the three body decay is dimensionless, Eq.(f47l) implies that the function 
F{E\,E 2 ) is also dimensionless and is related to the function F4 of |J by F(E%, E 2 ) = m£> s F±{s\, s 2 , Q 2 ■ 

m D s )- 

Combining Eqs.(|47|) and (|48|) , we get 

< Tr + ir~ir + \T\D+ > WA = % ai ° F ™ Ds V c * s V ud F(E 1 ,E 2 ) (49) 

V2 m Da 

The corresponding branching ratio is written in the form 

B{Dt -» tt + it-it + ) wa = N x / (50) 

where N is a normalization factor, 
and / a phase space integral given by 



I = —q— I I dE x dE 2 \F(E U E 2 )\ 2 (52) 



In Eq.(51) the factor i is again due to the presence of the two tt + in the final state. 



The value of N is computed using the same parameters as in section III and the result is N = 
1.0399 • 10- 2 . 

2). In the non resonant case we assume the structure function F to be a constant independent 
of the pion energies such as F(Ei, E 2 ) = Fnr. The experimental non resonant rate quoted in Eq.(|8|) 
allows us to extract the modulus of Fjy^. The corresponding phase space integral Inr in Eq.(|52|) with 
constant Fnr has the value 

I NR = 0.1053 \F NR \ 2 (53) 
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and we obtain by using Eqs.(g), (|50|) and (|53| ) 

\F NR \ = 3.04 i ° f s (54) 



3). A decomposition formally similar to Eq.(|49|) in terms of F(E\, E2) can be conveniently written 
for the two other components of Eq.([lFJ). This is obvious for the pir component which is of the W 
annihilation type and we get from Eqs.([T3D and (pl|) 



with 



Ffm {E l ,E 2 ) = A p7T ■ H pn (E X ,E 2 ) (55) 



Afm ~ 1 m 2 , -m 2 -im,T ,{m 2 ) 9 *'<™ 9p ™ ~ * % ^ (56) 

H^Ek) = 2 S2 ~ S ' r + 2 Sl ~ S " r (57) 
- si - « y/sl T p m z p - s 2 -i y/s^ T p 



For the fair component which is of the spectator type we obtain, from Eqs. (|3q) and (|39|), the structure 
function Ff o7r (Ei, E 2 ) rewritten in terms of the form factor F^ s ^° (m 2 ) 

Ff Tr(Ex, E 2 ) = A fon Hf Q7T (Ei, E 2 ) (58) 

with 

A '°' = fe) (&) ('"J) F » D - ,0 ( m -)*-- = i^-i < 59 > 

2 2 
H f0W (Ei,E 2 ) = — /0 . ^ r + — /0 . „ (60) 

In the spectator model the form factor is real, such that 4>f 07T = or tt. 
The total structure function F(Ei,E 2 ) has three components 

F(E U E 2 ) = F piT (Ei, E 2 ) + Ff 07T (Ei, E 2 ) + Fnr (61) 

and the computation of the total rate is made using the formula 

B{D+ -> 7r + 7T-7r + ) = 1.0399 • 10" 2 / (62) 



where I is the phase space integral defined in Eqs.(p2[). 

Let us first consider the phase space integrals involving the functions H{E\, E 2 ). We have previ- 
ously obtained 

' ' ' dE x dE 2 = 0.1053 (63) 



in 



1 



Hf o7T (Ei,E 2 )\ dEi dE 2 = K = 5.7217 (64) 
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m D a 



\H p7T (E 1 , E 2 )\ 2 dE 1 dE 2 = J = 8.0997 (65) 



For the interferences between the three components we need three more integrals 

H fon = -4- / / H foir (E 1 ,E 2 )dE 1 dE 2 = \H fQV \ e l<t>H ^ (66) 



m 



D. 



m D s 



Hpo^ = — H p7T (E 1 ,E 2 )dE l dE 2 = \H pn \ e* ^ (67) 



H pf = A- I I H; 7T (E ll E 2 )H f07T (E 1 ,E 2 ) dE 1 dE 2 = \H pf \ e l +*n (68) 



m D a 

The numerical results are 



H foA 


= 0.2677 


(t>Hf TT = 


89.6° 


(69) 


1 HpIT | 


= 5.58 • 10~ 5 


<t>H P TT = 


-84.5° 


(70) 


\ H pf\ 


= 0.3868 


<pHpf = 


241.5° 


(71) 



The numerical values of |-F/v_r| and |^4/ 07r | are taken from the measured rates respectively in Eqs.( 
and (||). In the pn case we write A pw = A A™® 00 with < A < 1 where A = 1 corresponds to the 
C. L. experimental limit in Eq.(|lO|) . 

There is a priori no reason for the non resonant amplitude Fnr to be real and we allow an arbitrary 
phase (f>NR- The two other phases cpf 07T and 4> p7T are determined by our model as 

(ppn = (70 ±10)° for rV = (400±200)MeF (72) 

hon = 0° for F^ fo (ml) > 0. (73) 

The formula for I(tt + tt^tt + )tot is then written as 



I(tT + TT TT + )tOT = /(7T + 7T 7T + )nR + 7T + ) J ()7r + A 2 I(lT + lT lt A 

+ 2 \F NR \ \A folT \ \H fow \ Cos[(pf 07r + 4>Hf 7r - 4>nr] 

+ 2 A \A pn \ max \A fon \ \H pf \ Cos[^ fo7T - ^ + <t> HpK ] (74) 

+ 2 A | F NR \\Apnl max | H p7T | Cos [0 pn + (f> Hfm - <j) NR ] 

The errors for /(7r + 7r _ 7r + )TOT are computed in quadrature from the experimental errors on /(7r + 7r _ 7r + ) 
and I(7r + 7T _ 7r + )j 07r and for the errors of the amplitudes \F^r\ and |^4/ 07r |, we use the approximate 
relation Ax ~ ^Ax 2 . 

The result of our calculation of I(tt + tt^tt + )tot directly gives the total branching ratio by the 
relation (^2|) and is represented in Figs. 3 as a function of the phase parameter 4>nr for the two 
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extreme cases A = and A = 1. The one standard deviation domains can be compared with the 
experimental result of Eq.(R). 

It is clear from Fig.3-a corresponding to a two component model (A = 0), and Fig.3-b where a max- 
imal contribution has been added (A = 1), that there is no difficulty in fitting the experimental 
result (H) within one standard deviation, for any value of the phase parameter 4>nr- 

VI. Energy distributions of the pion in — > 7r + 7r + 7r decay 

1) . The two 7r + meson energies being E\ and E2 with our convention of section II, the function 
F(E\,E 2 ) introduced in section V is symmetrical in the exchange of E\ and E 2 . Such a property 
extends to the Dalitz plot which is a double differential distribution proportional to \F{E\, -E^)! 2 - 
However it is probably premature to discuss the detailed properties of the Dalitz plot and it is more 
realistic to consider single differential quantities like the 7r meson energy distributions. 

We must keep in mind that a distribution in the pion energy Ej is just the mirror of a distribution 
in the invariant two pion mass squared Sj because of the relations ([!]), Sj = m 2 D + m 2 — 2rriD s Ej. 

2) . For the final state 7r + 7r + 7r~, we define the ir + and ir~ energy distributions in the following way 

iriDt = *°+< £ '> • dT(Dt 1e 3 + * + *~ ] = * a -^ < 75 > 

where the distributions G+(E\) and G-(E 3 ) are defined by 

G+CEi) = / \F(E 1 ,E 2 )\ 2 dE 2 (76) 

m Ds JE- (Ei) 

G-(E 3 ) = / \F(m Ds -E 2 -E 3 ,E 2 )\ 2 dE 2 (77) 

m Ds Je-(e 3 ) 

and the limits of integration E±(E) have been given in Eq.(||) of section II. 
The normalization constants N and M are simply related by 

M = ^f - m Ds M (78) 

and the numerical value of J\f is M = 7.446 • 1CP 15 . 

The 7r + energy distribution G+(Ei) and the tt~ energy distribution G-(E 3 ) for all 7r + 7r + 7r _ events 
have been computed performing a single variable integration with the relevant functions defined in 
section V, |#/ 07r (£i, E 2 )\ 2 , {H^Eu E 2 )\ 2 , H fon {E 1 ,E 2 ), H /m (E 1 ,E 2 ), and H* T (Ei, E 2 )Hf 07T (Ei, E 2 ). 
We observe that the last three integrals are complex with an energy dependent phase and, as a con- 
sequence, the nature of constructive or destructive of the various interferences is an energy dependent 
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concept. The terms Fnr, Af on , and A pn are treated as in section V with their moduli extracted from 
the measured branching ratios. The two extreme cases A = and A = 1 have been considered. 

It is clear that the shape of the tt meson energy distributions G+(E\) and G^{E^) depends on 
the constant phase of the non resonant amplitude F^r in a very sensitive way. We notice in Q 
that we have a Dalitz plot in the (si,S2) plane corresponding to the (68.1 ± 12.4) events observed in 

— > 7r + 7r~7r + decays. A projection of this Dalitz plot on the s± axis shows clearly a peak associated 
to the /ott + quasi two body state and an important background. Of course, such an histogram in Q 
is nothing but our ir + meson energy distribution G+(E\). 

The G+(Ei) distribution has been studied for values of 4>nr between 0° and 360°. Comparing our 
theoretical curves with the shape of the experimental histogram [|J , we obtain an acceptable agreement 
if the angle <Pnr is in the range 

70° < 4>nr < 140° (79) 

Let us now come back to Eq.(^9|) where the function F(E\,E2) is defined. Beside the phase of 
the CKM factors, we have the phase n/2 due to the factor i entering in the definition of the matrix 
element < |./L|-D+ > and the phase of the function F(Ei, E2). In the non resonant component, the 
last phase is <Pnr- We observe that the value 4>nr = 90° compatible with the range ( |79| ) makes the 
quantity i Fnr real and negative. In the absence of more detailed experimental informations we shall 
make a simple choice 4>nr = 90° in what follows. 

We have represented in Figs. 4 and Figs. 5 the pion energy distributions G+(E\) and G-(E%) for 
the two cases A = and A = 1 using (/)nr = 90°. For completeness we give the theoretical branching 
ratios corresponding to <Pnr = 90° which can be directly extracted from Figs. 3. 

(i) A = 0, B(D+ -» ir + ir-ir + ) TH = 2.022 ± 0.528 %, (80) 

(ii) A = 1, B(£>+ -» ir + Tr-ir + ) TH = 2.259 ± 0.524 %. (81) 

Both values Eqs.© and © agree with the experimental one of Eq.(|), (1.54 ±0.42) %. 

VII. Predictions for D+ — > 7r°7r°7r + decay mode. 

We have no experimental data on the decay mode Df — > 7r 7r 7r + . However we shall use the same 
three component model for this decay mode and let us start our discussion with isospin considerations 
concerning these three components. 

1). Consider first the case of a non resonant 7r°7r°7r + state. The constant Fnr has no reason 
to be the same for the two modes tt + it~it + and 7r°7r°7r + . By assumption, with a constant function 
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Fnr(Ei, E2), we have a full symmetry in space between the three pions. From the Bose-Einstein 
symmetry, the isospin configuration has also to be totally symmetric. Consider now a third rank fully 
symmetric tensor in a three dimensional space. It has 10 independent components. With respect 
to the isospin SO (3) orthogonal group, such a tensor is reducible into an isospin 1 = 3 part with 7 
components and an isospin 1 = 1 part with 3 components. In our specific W annihilation model for 
the non resonant three pions, only the later part contributes, the ud weak current being an isovector. 
By inspection of the relevant Clebsch-Gordan coefficients, we get the result 

^(vr+Tr+vr-) = 3 F nr (tt\\ + ) . (82) 

As a consequence, in our model, we obtain for the non resonant part : 

rp+ _> nW^R 1 11-1 0/ m) 

JW^WW " 9 - 1L1 % ' (83) 

and the non resonant branching ratio B(Df — > 7r°7r°7r + )A^ is expected to occur only at the 10~ 3 level. 

2) . The quasi two body state p°ir + cannot produce a 7r°7r°7r + final state because of the isovector 
character of the p meson. However p + 7r° can give such a state when followed by the decay p + — > 7r + 7r°. 
The Df meson being an isoscalar, the decay amplitudes for D s — > p°ir + and D s — > p + 7r° are equal. 
When mass differences between charged and neutral p's and 7r's are neglected we obtain 

T(Df^p\ + ) = T(Df^p + 7r°) (84) 

With equal decay widths for p° — > 7r + 7r~ and p + — > 7r + 7r° we obtain the ratio of rates 

r(D+ p%+ ^+^-^+) 1 j 

We observe that the finite p width corrections are the same in both final states 7r + 7r~7r + and 7r°7r°7r + 
because two Breit Wigner amplitudes are present in the variables si and §2 with the choice made in 
section II. 

3) . For the third component, i.e., the decay — > /qtt + followed by /o — > 7r + 7r _ (/o — ► 7r°7r°) 
produces a 7r + 7r _ 7r + (7r 7r 7r + ) final state. Because of the isoscalar character of the /o(980) we have 

r(/ -> 7T+7T-) = 2 r(/ ^ /V ) (86) 

On the other hand, in the Df — > /o7r + — > 7r + 7r _ 7r + case, we have two Breit- Wigner amplitudes in 
the variables Si and S2 whereas in the — > /ovr + — > 7r°7r°7r + case there exists only one Breit- 
Wigner amplitude in the variable S3 with the choice made in section II. As a consequence the phase 
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space integral K introduced in Eq.(|4l|) has the value K c = 5.7217 for the final state ir + ir tt + and 
Kpf = 2.7469 for the final state 7r 7r°7r + , and we obtain the ratio of rates 

r( J P+^/ovr+^vr A+) 1 K N 

r(£>+ - / o7r + - 7r+7r-7r+) " 2 K c " 24/ ° 

The departure of this ratio from j is simply due to the interference between the two Breit-Wigner 
contributions in the Df — ► /ovr + — > 7r + 7r~7r + case. 

4) . The function F(E\,E 2 ) associated to the final state 7r°7r°7r + is written in the form 

F(Ei,E 2 ) = -F NR + —j=Af Q1T Hf o7r (Ex, E 2 ) + \ A™ x H pn (Ei,E 2 ) (88) 

where in the right hand side of Eq.(|8S|) all terms but H t 07r (Ei, E 2 ) are the same as those defined in 
section V for the tt + it~tt + final state. In Eq. (|88|) , E\ and E 2 are the 7r° energies and £"3 the ir + energy. 
Whereas Hf on (Ei, E 2 ) in the 7r + 7r~7r + case contains two Breit-Wigner terms with the variables s\ and 

s 2 (Eq.(|60|) ), for the 7r°7r°7r + final state we have only one Breit-Wigner term with the variable S3 

2 
fa 

Hf on (E 1 ,E 2 ) = — g : — ; — ^ ( 89 ) 

The total rate is now computed by integration the quantity \F{E\,E 2 )\ over E\ and E 2 . The 
procedure of calculation is similar to the one explained in detail in section V and the result for 
I ^7r°7r°7r + ) tot is shown on Figs. 6 for the two extreme cases A = and A = 1. Because of the isospin 
factor 3 in front of Fnr the role played by the non resonant component is obviously less important 
in the 7r°7r°7r + case than it was in the 7r + 7r~7r + one. For such a reason the quantity I(7r 7r 7r + )tot m 
the two component model, A = 0, is insensitive to the phase parameter 4>nr- Again due to the isospin 
factors, the possible role of the quasi two body state pn might be more important in the 7r°7r°7r + case 
than it was for the tt + tt^tt + one. Such a qualitative expectation is clearly seen in Figs. 6. Both Figs. 3 
and 6 are drawn at the same scale and the order of magnitude of the branching ratio B{Df — > 7r°7r°7r+) 
is expected to be few 10~ 3 and that might be the reason why the mode Df — > 7r°7r°7r + has not yet 
been experimentally observed. 

With the choice 4> NR = 90° we find 

(*) A = 0, B{Dt -> ^°vr 7r + ) T ^ = 0.294 ± 0.092 %, (90) 

(m) A = 1, B{Dt -> -K -K°it + ) T H = 0.604 ± 0.101 %. (91) 

5) . The vr° and 7r + energy distributions Gq{E\) and G+(E$) for the final state 7r°7r°7r + are 
computed with the formulae (|76|) and (f77| ) using now the function F(E\,E 2 ) given in Eq.(^). The 
procedure of calculation is similar to the one explained in section VI for the Tr + ir~ir + final state. 
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We have represented respectively in Figs. 7 and 8 the ir° and ir + energy distributions Gq(Ei) and 
G+(Ez) for the two extreme situations A = and A = 1, using as previously the phase cpNR = 90° for 
the non resonant amplitude. Figs. 8 show clearly the peak due to the quasi two body fo7r + state and 
a possible quasi two body p + 7r° contribution is clearly seen in Fig.7-b. 

Let us end by considering the tt meson energy distribution corresponding to an uniformly populated 
phase space of a constant amplitude. The energy distributions for each pion are obviously identical 
and in the case — > (tt + tt^tt + )n R we obtain, using Eq.([|), the analytic form 

, f (E 2 - ml)(m 2 - 3m 2 - 2m D E) 1 1/2 
G NR (E) = \F NR \ 2 \ K - i\ Ds 2 J-, \ 92 

Using the experimental value of \Fnr\ 2 , we obtain the distribution represented in Fig. 9. 



VIII. Summary and Concluding Remarks 

We have studied, in this paper, the decay of the Dj~ meson into three n mesons by analyzing the 
experimental data available for the decay mode Df — > ir + Tr~ir + . The amplitude can be considered as 
the superposition of two mechanisms : the spectator decay of the c quark contained in the Df meson 
and the W + annililation into three pions. Both mechanisms are needed to explain the data and the 
latter appears to be considerably important. 

The non observation of the pir quasi two body state allows us to determine an upper limit of about 
10 MeV for the leptonic decay constant / / of the it' = 7r(1300) meson. Such an upper bound of / / 
is compatible with theoretical estimates using non-relativistic chiral quark model |l(J and QCD sum 
rule technique |11]. We remark that the commonly used values of the f > decay constant [[?| are much 



larger than our upper limit. 

The large /ovr + branching ratio compared to the p°ir + one allows us to exclude the dominant 
(uu + dd)/y/2 quark structure for the /o(980) meson. Assuming the /o(980) to be a ss state, the 
spectator mechanism will govern the decay — > /o(980)7r + and from experimental data it is possible 
to estimate the hadronic form factor F^ B ^°(m 2 ) to be 0.36 t o^g- 

Furthermore, a large branching ratio of the meson into a non resonant three pion state has 
been experimentally measured. This decay can only proceed through the WA and the corresponding 
decay amplitude is assumed to be independent of the n meson energies. We then introduce a phe- 
nomenological complex constant -F/vr, its modulus is extracted from the non resonant branching ratio 
B(D+ -> 7 r+7r-7r+) Jvfl . 

For the full decay amplitude, we propose a three component model involving the non resonant 
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background, the fo^ + quasi two body state, and a possible pit quasi two body state with a parameter 
A such that for A = 1, the branching ratio B(Df — * 

A + ) 

takes the value of the 90% confidence level 
upper limit. We first compute the total rate for Df —* 7r + 7r~7r + as a function of the phase 4>nr of the 
constant amplitude Fnr. Agreement with experiment is obtained, within one standard deviation, in 
both cases A = and A = 1, for all values of (j>NR between 0° and 360°. 

For the decay mode Df — * 7r + 7r~7r + , we compute the tt + and tt~ energy distributions which 
depend on the phase parameter (ft nr.- A qualitative agreement between our ir + distribution with the 
experimental histogram in the invariant (ir + ir~) mass is obtained when 4>nr is restricted in the 
range given by Eq.(fF^) including 4>nr = 90°. Fixing now 4>nr = 90°, we present our predictions for 
the 7r + and tt~ energy distributions including one standard deviation errors deduced for those of the 
experimental rates. 

The same model is used to make preditions for the decay mode Df — ► 7r°7r°7r + which is not 
yet experimentally observed. The total rate for Df — ► 7r°7r°7r + is found to be essentially one order 
of magnitude smaller than the Df — * ir + ir~ir + rate and this might be the reason why the decay 
Df — > 7r°7r°7r + has not yet been detected. Also the 7r° and tt + energy distributions are predicted by 
our model. 

If the quasi two body states fair and pir can be theoretically understood, the non resonant part 
has been treated, in this paper, only from a purely phenomenological view point, using a constant 
complex parameter Fnr- It would be interesting to explain theoretically the experimental values 
found for |-F/vr| and 4>nr- To our knowledge no model has been proposed for such case. 

The detailed understanding of the WA is very important in its own right, and for that purpose, 
the Df, like the Bf , are particularly auspicious for WA to manifest in its full strength due to both 
color and CKM favored factors. In all other cases, i.e., D+,D°,B + ,B°, and B°, the WA is either 
masked by the color and CKM suppressed factors, or contaminated by the spectator decay mechanism 
at the quark level, such that unambiguous WA effect can be hardly isolated. Therefore experimental 
and theoretical investigations of the Df ( and the Bf later) decays into pions, in both inclusive 
and exclusive modes, are of great interest. They are the ideal laboratories for studying the weak 
annihilation mechanism. 
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Figure captions 



1. Figure 1 : W annihilation mechanism for the Df — > p° + tt + decay mode. 

2. Figure 2 : Spectator mechanism for the Df — > /o + n + decay mode. 

3. Figure 3 : Total branching ratio for the decay Df — > tt + tt^tt + as a function of the phase 
parameter 4*nr, (a) hi the two component model : A = 0, (b) in the three component model 
with the maximal pn + contribution : A = 1. One standard deviation errors are indicated and 
the horizontal bound is the experimental result including one standard deviation errors. 

4. Figure 4 : The ir + meson energy distribution function G + (E\) for the hnal state ir + n~n + , 
(a) for A = in the two component model, (b) for A = 1 in the three component model with 
the maximal pir + contribution. One standard deviation errors are indicated. Here the quantity 
Ef = 0.7453 GeV is associated to the /o resonance and E p = 0.8386 GeV is associated to the 
p° resonance. 

5. Figure 5 : The 7r~ meson energy distribution function G-(Es) for the final state it + it~it + , (a) 
for A = in the two component model, (b) for A = 1 in the three component model with the 
maximal pir + contribution. One standard deviation errors are indicated. 

6. Figure 6 : Total branching ratio for the decay Df — ► 7r°7r°7r + as a function of the phase 
parameter 4>nr, (a) hi the two component model : A = 0, (b) in the three component model 
with the maximal pir + contribution : A = 1. 

7. Figure 7 : The it meson energy distribution function Gq(Ei) for the final state 7r°7r°7r + , (a) 
for A = in the two component model, (b) for A = 1 in the three component model with the 
maximal p + ir° contribution. One standard deviation errors are indicated. Here the quantity 
E p = 0.8386 GeV is associated to the p resonance. 

8. Figure 8 : The tt + meson energy distribution function G + {E^) for the final state 7r°7r°7r + , (a) 
for A = in the two component model, (b) for A = 1 in the three component model with the 
maximal p + 7r° contribution. One standard deviation errors are indicated. Here the quantity 
Ef = 0.7453 GeV is associated to the /o resonance. 

9. Figure 9 : The ir meson energy distribution Gnr(E) for the non resonant part of the decay 
Df — ► 7r + 7r _ 7r + . One standard deviation errors are indicated. 
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Figure 1: W annihilation mechanism for Df — >■ p + 7r + decay mo 
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Figure 2: Spectator mechanism for Df — > / + 7r + decay mode. 
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Figure 3: Total branching ratio for the decay D+ — > ■k + -k~-k + as a function of the phase parameter 
4>nr, (a) in the two component model : A = 0, (b) in the three component model with the maximal 
pn + contribution : A = 1. One standard deviation errors are indicated and the horizontal bound is 
the experimental result including one standard deviation errors. 
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Figure 4: The n + meson energy distribution function G + (E{) for the final state 7r + 7r _ 7r + , (a) for 
A = in the two component model, (b) for A = 1 in the three component model with the maximal 
pn + contribution. One standard deviation errors are indicated. Here the quantity E fo = 0.7453 GeV 
is associated to the f resonance and E p = 0.8386 GeV is associated to the p° resonance. 




Figure 5: The n~ meson energy distribution function G_(E 3 ) for the final state 7r + 7r~7r + , (a) for 
A = in the two component model, (b) for A = 1 in the three component model with the maximal 
pn + contribution. One standard deviation errors are indicated. 
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Figure 1: Total branching ratio for the decay Df — > 7r°7r°7r + as a function of the phase parameter 
4>nr, (a) in the two component model : A = 0, (b) in the three component model with the maximal 
pn + contribution : A = 1. 
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Figure 2: The 7r° meson energy distribution function G (E 1 ) for the final state 7r°7r°7r + , (a) for A = 
in the two component model, (b) for A = 1 in the three component model with the maximal p + 7r° 
contribution. One standard deviation errors are indicated. Here the quantity E p = 0.8386 GeV is 
associated to the p resonance. 
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Figure 3: The n + meson energy distribution function G + (E 3 ) for the final state 7r°7r°7r + , (a) for A = 
in the two component model, (b) for A = 1 in the three component model with the maximal p + 7r° 
contribution. One standard deviation errors are indicated. Here the quantity E fa = 0.7453 GeV is 
associated to the f resonance. 
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Figure 4: The n meson energy distribution G NR (E) for the non resonant part of the decay D+ — > 
7T + 7r~7r + . One standard deviation errors are indicated. 



